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A new generation of angular-resolved photoemission spectroscopy (ARPES) measurements on the 
cuprate superconductors offer the promise of enhanced momentum and energy resolution. In partic- 
ular, the energy and temperature dependence of the on-shell nodal (k x — k y ) quasiparticle scattering 
rate can be studied. In the superconducting state, low temperature transport measurements suggest 
that one can describe nodal quasiparticles within the framework of a BCS d-wave model by includ- 
ing forward elastic scattering and spin-fluctuation inelastic scattering. Here, using this model, we 
calculate the temperature and frequency dependence of the on-shell nodal quasiparticle scattering 
rate in the superconducting state which determines the momentum width of the ARPES momentum 
distribution curves. For a zero-energy quasiparticle at the nodal momentum fcjv, both the elastic 
and inelastic scattering rate show a sudden decrease as the temperature drops below T c , reflecting 
the onset of the gap amplitude. At low temperatures the scattering rate decreases as T 3 and ap- 
proaches a zero temperature value determined by the elastic impurity scattering. For T > T c , we 
find a quasilinear dependence on T. At low reduced temperatures, the elastic scattering rate for 
the nodal quasiparticles exhibits a quasilinear increase at low energy uj which arises from elastic 
scattering processes. The inelastic spin-fluctuation scattering leads to a low energy uj 3 dependence 
which, for uj > 3Ao, crosses over to a quasilinear behavior. 

PACS numbers: 74.25.Jb, 74.20.Fg, 74.72.-h 
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I. INTRODUCTION 

A d-wave BCS framework has proved useful in describ- 
ing the thermal and microwave Q conductivities of 
the superconducting cuprates 0,11, El- This suggests that 
a quasiparticle description of the nodal excitations is ad- 
equate, at least at low energies. Here, the frequencies are 
small compared with the gap magnitude and the temper- 
ature dependence of the nodal quasiparticle lifetime plays 
the dominant role. Recent angle-resolved photoemission 
spectroscopy (ARPES) experiments on Bi 2 Sr 2 CaiCu208 
(BSCCO) are providing high resolution data on the mo- 
mentum and energy dependence of the nodal quasipar- 
ticle lifetimes as well as their temperature dependence 
in this low energy range H, 0, El- Thus, the ques- 
tion arises whether transport and ARPES data can be 
understood within a single framework. Establishing this 
agreement is important to confirming the nature of the 
superconducting phase as a BCS-like d-wave state with 
nodal quasiparticle excitations. 

The width Ak(uj,T) of the momentum distribution 
curve (MDC) measured by ARPES experiments is pro- 
portional to the inverse of the on-shell quasiparticle life- 
time. Along a nodal cut in momentum, with k x = k y , 

the magnitude k = \/k x + k'l is set by ui = et, i.e. k 
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becomes effectively a function of uj; we refer to this as 
"on-shell" . Here, Sk is the quasiparticle energy. Thus, 
the lifetime of such a nodal quasiparticle depends upon 
uj and the temperature T. As k decreases from the nodal 
Fermi momentum kjq, uj can be swept over an energy 
range which is significantly larger than Ao. Therefore, 
there is a need to extend the inelastic scattering lifetime 
calculations to cover a wider range of energies. In addi- 
tion, impurity scattering, particularly due to out-of-plane 
forward scattering, is believed to be important so that it 
is also of interest to determine its temperature and en- 
ergy dependence 0, 0] . 

In the following, we will find it convenient to discuss 
the elastic and inelastic quasiparticle scattering in terms 
of a scattering rate 



r(w,k,T) = -£"(w,k,T) 



(1) 



where S" is the imaginary part of the quasiparticle 
self-energy. The inverse of the quasiparticle lifetime 
l/r(uj,k,T) = 2r(w,k,T) and the width of the MDC 
about a given k determined from Eu = uj is 



Ak(uj,T) 



2r(q,,kU= ek ,T) 

«(k) 



(2) 



with v(k) = the bare band velocity. Here, we re- 
port on results obtained within a d-wave BCS framework 
for the energy and temperature dependence of the nodal 
quasiparticle elastic and inelastic scattering rates and dis- 
cuss the uj and T dependence of T(u>, k| w=£fc , T). 

We begin in Section II with a discussion of the temper- 
ature and energy dependence of the elastic impurity scat- 
tering contribution. In BSCCO, there are both in-plane 
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and out-of-plane scattering centers The in-plane 

impurities are believed to give rise to strong scattering 
and are often treated in the unitary limit [lj]. On the 
other hand, the out-of-plane impurity scattering is weak 
and tends to be forward. It will be treated within a self- 
consistent Born approximation. We will see that the tem- 
perature dependence of the elastic scattering rate arises 
from the opening of the gap as T decreases below T c , 
while its energy dependence reflects the decreased phase 
space as u> becomes smaller than the gap magnitude. The 
low energy inelastic scattering discussed in Section III is 
dominated by short range Coulomb scattering. At low 
energies, the linear uj dependence of the density of states 
gives rise to an inelastic scattering rate which varies as T 3 
for T larger than ui or u> 3 for lu larger than T. At higher 
energies, one can probe the w-depcndcncc of the effective 
interaction. Here, we confine ourselves to the Hubbard 
model and explore the spin-fluctuation contributions to 
the higher energy behavior of the inelastic scattering. In 
Section IV, we discuss the combined effects of the elastic 
and inelastic scattering processes, and Section V contains 
our conclusions. 



II. ELASTIC SCATTERING 

Elastic scattering in BSCCO can arise from both im- 
purities and disorder in the CuO planes as well as from 
regions outside these planes. The scattering rate due to 
unitary scatterers (possibly Cu vacancies) with concen- 
tration roughly n u ~ 0.2% observed as zero-bias reso- 
nances in scanning tunneling microscopy (STM) exper- 
iments i s we H understood theoretically. It will be 
treated as usual in the self-consistent T-matrix approxi- 
mation, 



-eLu 



E k G(k, 



TO, 



(3) 



where tj are the 2x2 Pauli matrices in Nambu nota- 
tion. At low frequencies and temperatures, the unitary 
scatterers give rise to an impurity resonance of height 
Klu(^ = 0) = lu v/T^Ao" ~ 10" 2 t, and roughly the 
same width, with t the nearest-neighbor hopping matrix 
element. In the normal state we find a scattering rate 
of T u ~ n u Ep ~ 10 _3 i due to these defects, leading to 
7 U ~ 1 to 2 meV. If, as we believe, the normal state elas- 
tic scattering rate at the node is a significant fraction of 
the gap, the contribution of the unitary scatterers will be 
difficult to resolve in ARPES experiments on BSCCO. 

The elastic scattering from out-of-plane impurities, 
such as interstitial O ions, can be modeled by a 
momentum-dependent potential, 



V(k,k') =V f(k,k') 



(4) 



which has a forward scattering form factor /(k, k') which 
cuts off the scattering when |k— k'| exceeds a momentum 
K, as has been discussed in Rcf. [EJ We will measure k 



in units of kp so that k _1 characterizes the range of the 
impurity potential in units proportional to the Cu-Cu 
spacing a. Since the scattering from an out-of-plane im- 
purity is relatively weak, it can be treated within the self- 
consistent Born approximation [l5j . For a nodal quasi- 
particle with k x = k y , the gap vanishes and the scat- 
tering rate due to elastic scattering is determined from 
the imaginary part of the To Nambu self-energy in the 
self-consistent Born approximation, 



E ei , (k, W ) = r il V|^(k,k 



/\|2 



k' 



(5) 



where is the planar density of the out-of-plane im- 
purities. In addition, because the potential V(k, k') is 
anisotropic, it is necessary to renormalize the order pa- 
rameter as well: 



E rfl i(k,w)=n i 2JV(k,k')| 
fe' 



A k - 



f 2 - A 2 



(0) 



where a) =u> — E^o and = Ak + E e ;.i- The renormal- 
ization of £k by disorder is also nonzero in general but 
was shown to be negligible in Ref. Il2l and will therefore 
be neglected here. For our numerical calculations we will 
use a simple parameterization of the band structure with 



e fc = — 2t (cos k x + cos k y ) — At' cos k x cos k y — fi (7) 



and 



A fc = ^ (cos k x - cos k y ) . 



(8) 



Here, t' = —0.35 and fi = —1.1, with energy measured in 
units of nearest neighbor hopping t. For the temperature 
dependence of the gap we use a common interpolation 
form 



(9) 



with a = 3, 2A / kT c = 6, and A = 0.2. These pa- 
rameters are chosen to give the observed magnitude of 
the T = gap and a somewhat more rapid rise at T c , 
consistent with experiment. 

The scattering rate r e ;(kjv,o;) is now determined by 
the imaginary part of the solution of Eq. (0 at the nodal 
wavevector k^r. Before turning to a numerical evaluation 
of Eq. ©-I©, we consider a simple analytic approxima- 
tion similar in spirit to that discussed in |l7j | in the con- 
text of T c suppression by disorder Measuring the 
momentum transfer from k^v along the k± and fen coor- 
dinates shown in Fig. ^ the imaginary part of Eq. (JSJ 
can be written as 



T el = -m\V \ 2 Im 



dk'dk' 



(2n) 



£ l/(k,k') 



v\ k'l 



v 2 k' 2 



(10) 
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FIG. 1: (Color online) The Fermi surface corresponding to 
the band Sk given by Eq. Q with t'/t = —0.3 and /i = —1.1. 
The difference of the scattered wave vector k' from the nodal 
wave vector kjv is shown as its parallel displacement along 
the Fermi surface kii and the perpendicular displacement k±. 
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FIG. 2: Gap slope renormalization V2/V2 vs. inverse scattering 
range k from the solution of Eqs. il l 21 and 1)14^ . 



Here, v± = vp at the nodal point and v 2 is the renormal- 
ized gap velocity. We have assumed that the k cut-off 
prevents scattering to other nodes. 

Since v\ 3> i>2, the n cut-off primarily affects the k\\ 
integral and one can integrate freely over the important 
range k± giving 



r ei = ^ ni N \V \ 2 Re 



dk'u 



vlkf 



(11) 



with Nq = {2ttvf)~ 1 the band density of states in units 
in which = 1. Setting Co = u> + iT e i, i.e. neglecting 
the real part of the self-energy, this becomes 



Rc 



V 2 K 



(uj + iT ei ) sin 1 



v 2 k 



i T P 



(12) 



where To = m Nq |Vo| 2 k is the normal state elastic scat- 
tering rate and v 2 k is the maximum renormalized order 
parameter probed by the scattering. From this result 
we recover immediately the uj — > 0, ku — > limit that 
T e i -> T from Ref. d 

The gap velocity at the node in the pure system, v 2 — 
2Ao/fcw, is significantly renormalized in the presence of 
disorder to v 2 . From Eq. ijfJJ we find, for fcii k, 



(v 2 - v 2 )k\\ 



leading to 



£0 

2k 



Im 



dk' H 



U2(A[| + fc||) 



as) 



v 2 - v 2 



v 2 



(14) 



This result, which should be solved self-consistently to- 
gether with Eq. I|12|) , obtains only in the limit where the 
scattering is sufficiently forward such that the momentum 
integration takes place only near the node; it interpolates 
correctly to the isotropic limit k — > 00, however, where 



the gap velocity renormalization vanishes, v 2 = v 2 . In 
the forward scattering limit k — ► 0, it is easy to check 
that v 2 — > v 2 /2, and the dependence for small k is oth- 
erwise weak, as shown in Fig. |21 

We will consider two different ways of characterizing 
the scattering rates of quasiparticles with k x = k y . First, 
suppose k is fixed exactly at the nodal Fermi surface 
kj\r. In this case the quasiparticle energy to = and 
one can study T e i(k.N,uj = 0) as a function of temper- 
ature. The contribution of the out-of-plane elastic im- 
purity scattering to the temperature dependence of the 
MDC line width Ak(T) for a quasiparticle with k = kyv 
is given by 2T el (k N ,uj = 0)/v F . When w = 0, Eq. 
can be written as 



1 = 



V 2 K 



hi 



1 + y/1 + (T el /v 2 ny 

(Tel/V 2 K) 



(15) 



At low temperatures, if v 2 k/Tq 1, we find the familiar 
Born result 



r e i(T) ~i> 2 Ke 



(16) 



with the gap maximum Ao in the isotropic case replaced 
by v 2 k. When v 2 k/Tq <C 1, 



(17) 



This limit, of course, always applies when T goes to T c , 
but it can also apply for all values of T/T c if the system is 
sufficiently dirty. Note that since To varies as |Vb| 2 K, if Vq 
is kept fixed, the normal state scattering rate decreases 
as k becomes smaller. This simply reflects the reduction 
of phase space as the scattering is restricted to be in a 
more narrow forward cone. At the same time, the drop 
in the scattering rate as the gap Ao(T) opens for T < T c , 
is smaller for reduced values of k. This reflects the fact 
that the maximum effective gap reached by a scattered 
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FIG. 3: Evaluation of the elastic nodal scattering rate T e i 
from Eqs. 1121 and 1141 . a) temperature dependence of 
r c i(ui = 0, T) for three values of the impurity potential range 
parameter k = 0.5 (solid), 1.0 (dashed), and 2.0 (dotted) with 
To = 0.25Ao; b) energy dependence of r o i(aj,T = 0) for the 
same parameters; c),d) same as a),b) but for F = 0.5A . 



nodal quasiparticle is reduced by k so that the opening 
of the gap is less effective in suppressing the scattering 
rate. Results obtained from solving Eqs. (|12|l and 114|) 
for To/Ao = 0.5 and 0.25 are shown in Fig.|3|a) and c) 
as a function of k. The temperature dependence of the 
clastic scattering arises from the temperature dependence 
of the d-wave gap, and therefore varies as T 3 at low T. 
Note that the interpolation form, Eq. @, actually gives 
a qualitatively incorrect result at low temperatures for 
a d-wave superconductor. Eq. @ implies an activated 
thermal depletion of the condensate, whereas the correct 
solution to the d-wave gap equation has a leading term 
varying as T 3 0, giving A (T) ~ A (l - bT 3 /A 3 ), 
where b is a constant of order unity. In the clean limit, 
it follows from Eq. I|16|) that 

r.,(T) * r.,(r = o) (i + ^ ~ T ' } t> 
"W°> ( 1 + 2 rS) 

while in the dirty limit, Eq. (|17(l implies that 

(v2k) 2 \ b(v 2 K) 2 



(18) 



r el (T) ~ r i 



r i 



erg 

6T 2 



3A 3 r 

bn 2 



3A rg 



rp6 



(19) 



In Eqs. (HHJ) and ^ we have assumed V2 — 2Aq and 
used the asymptotic results for V2/V2 derived above to 
obtain simple expressions for the clean and dirty limits, 
respectively. 

It is also interesting to study the lo energy dependence 
of T e i at a fixed temperature. Within the current ap- 



proximate framework represented by Eqs. I|12|) and 114|l . 
we will study r e i(w, kjy, T). For T small compared with 
T c , one obtains the results shown in Fig. El b) and d). 
The initial increase in the scattering rate as the bind- 
ing energy \lu\ increases reflects the opening of the nodal 
phase space with increased quasiparticle energy. When 
is large compared to V2K, the scattering saturates at 
its normal state value. By expanding Eq. (|12f> . we find 
that 



r e i(w) 




(20) 



so this saturation actually occurs from values of T e i some- 
what above To, as seen in Fig. El 

For forward scattering, where V2k/Tq is small, one can 
obtain a useful estimate of the quasilinear slope. In this 
case, r e i(w) rises from its value of r o i(o; = 0) at u> = 
to Tq over a frequency range of order v 2 k with a quasi- 
linear slope (ro — r e i(w = 0))/(v 2 k). In the dirty limit, 
Eq. {I7|), for < |cj| < v 2 k this gives 



r c ,M ~ r 0l (a, = o) + i 



(21) 



Finally, if the weak scattering is sufficiently isotropic 
(k > 1), we recover the standard result for pointlike weak 
scattcrers 

W ^r el (. = 0) + |^«^. (22) 

Which of the results, Eq. (|2T|) or (|2~2"|l is valid depends 
upon the character of the disorder. 

To check these estimates, and evaluate the scattering 
rate for various scattering potential ranges and impurity 
concentrations, as well as to treat momenta away from 
the Fermi surface more accurately, we introduce a more 
realistic model for the elastic scattering. We consider a 
screened exponential fall-off in 2D such that, 



|V(k,k')| : 



q 2 + k 2 



(23) 



where q = k k' is the momentum transfer. We then per- 
form a self-consistent solution of equations (JSJ and 10. 
Results for r c i analogous to Fig. [21 but taken "on-shell" 
at momentum k such that £k = w are presented in Fig. 01 
Results have been scaled such that the (self-consistent) 
normal state scattering rate is the same for all curves; 
as discussed above, this implies that the more forward 
scattering cases correspond to larger impurity concentra- 
tions. It is seen that the qualitative behavior close to the 
Fermi level is captured rather well by the approximate 
model discussed above where the scattering is restricted 
to a cone of width ~ k. It is also interesting to note 
that the quasilinear behavior at the low to values per- 
sists until extremely long scattering ranges, of order 5-10 
lattice spacings, for reasonable total normal state scat- 
tering rates. This suggests on the one hand that elastic 



5 



r ,(k M ,T )=0.5A n a) 

er N c' 



0.6 



C) 



u 



0.5 



0.3 



0.3 



O.li 



with 



0.5 
T/T 



r ,(kT )=A„ b) 
el v N c' 



0.9 



0.8- 



0.7„ L 



0.5 
T/T 




FIG. 4: Elastic scattering rate r e ] for the Yukawa-type po- 
tential with range k _1 Eq. I|230 . a): F c i/Ao vs. T/T c for a 
normal state elastic scattering rate r e i(fcjv,T c ) = 0.5Aq and 
k = 0.1, 1, and oo (isotropic), c): r c i/Ao vs. uj/Ao for the 
same k and r e i(fcjv,T c ) values, b) and d): same quantities 
plotted as above, but with r e i(fcjv,T c ) = Ao. Results for dif- 
ferent values of the scattering parameter « = 0.1 (solid), 1 
(dashed), and the isotropic case k = oo (dotted) are shown. 



scattering is a likely explanation for the quasilinear be- 
havior seen at low energies in the MDC width; on the 
other hand, it also means that it will be difficult to de- 
termine the range of the scatterers precisely from such a 
measurement alone. 



III. INELASTIC SCATTERING 

In this section we discuss the scattering rates due to 
inelastic electron-electron scattering by a phenomcno- 
logical short range Coulomb interaction based on the 
two dimensional Hubbard model. For the tight-binding 
bandstructure considered here, these scattering pro- 
cesses are dominated by exchange of antiferromagnetic 
spin-fluctuations. We are using a conventional Bcrk- 
Schrieffer-like theory as has been used before to 
discuss the single-particle inelastic lifetime [2(| and has 
proved to give a qualitative description of low-energy mi- 
crowave and thermal conductivity lifetimes |2ll as 
well as NMR relaxation rates |2^| in the cuprates. 

In second-order perturbation theory for the two dimen- 
sional Hubbard model, the imaginary part of the nodal 
quasiparticle self-energy due to inelastic scattering from 
the onsite Coulomb interaction U can be written as 



r incl (w,k,T) = -£f ne ;(w,k,T) = 



dfl [n(fi) + f(n - u)] 



xZ{q,n)N(k-q,u,-n) 



(24) 




(25) 



f(E k+q ) - f(E k ) 



Ek+qEk 



Ek+qEk 



— (Ek+q — Ek) " 

f(Ek + q) - f(E k ) 
Q + (Ek + q + E k ) + l0 + 

f(E k+q ) + f(Ek)-l 
Q-(Ek+q+E k ) + i0+ 



Here, Ek = + A 2 ,. As before, we are considering a 
nodal quasiparticle with k x — k y and the magnitude of 
k is set by u> — Ek- At low energies, where uj and T are 
small compared to the zero temperature gap amplitude, 
Eq. (|24|1 determines the form of the u> and T dependence 
of the inelastic quasiparticle scattering. Higher order pro- 
cesses lead to a renormalization of U, but the phase space 
restrictions imposed by the d-wave nature of the super- 
conducting state determine the low T and u> dependence 

of r inel (w,k| w=Ek ,T) damn. 

For k = kjy and u> = 0, the low temperature inelastic 
scattering varies as 



IWO.T) 



j>3 

A 5 



(27) 



with a prefactor of order one |20t l24j . For the set of 
parameters discussed below the prefactor becomes 2.4. 
It has been ar gue d that Fi ne i should vary as T 5 / 2 at low 
temperatures [25J . However, for the temperatures T > 
0.05T C we have studied, the T 3 behavior provides a much 
better fit to our numerical results in Fig. E{a). Similarly, 
at T = the to dependence is given by 



i(w,0) 



A 2 



(28) 



The extra power T 3 and uj 3 versus the T 2 and u 2 depen- 
dence of the usual (3D) Fermi liquid inelastic quasipar- 
ticle Coulomb scattering is due to the uj/A dependence 
of the d-wave density of states. Note that the usual T 3 
and ui 3 phonon contributions to the low energy quasi- 
particle scattering in a normal metal give contributions 
which vary as T 4 and u> 4 in a d-wave superconductor. 

At higher energies, the problem becomes more com- 
plicated. New, collective channels may open such as the 
so-called 7r-resonance or possibly the B\ g phonon. Here, 
we replace xo(<l,u) by the RPA form 



1 - Uxo(q,u) 



(29) 



and consider the spin-fluctuation contributions to the 
scattering of the nodal quasiparticle. Following the usual 
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FIG. 6: Left panel: Temperature dependence of inelas- 
tic RPA spin fluctuation scattering rate ri ne i(0, kjv, T)/Ao 
vs. T/T c for U — 2.2t. Right panel: Frequency depen- 
dence ri n ei(a;, k|^ =£k , T)/Ao vs. w/Ao for reduced temper- 
atures r/T c =0.1, 0.8, and 1.0. 



FIG. 5: Comparison of on-shell spin- fluctuation inelastic scat- 
tering rates for Yi ne i (uj, k|^ =£k , T)/Ao vs. lj/Aq for a reduced 
temperature T/T c = 0.1 . Dashed line: 2nd order perturba- 
tion theory with U = 6.7t. Solid line: RPA theory with 
U = 2.2t. Inset shows magnified low energy region. 



spin-fluctuation notation, we replace the coupling U 2 in 
Eq. HH) by §C/ 2 which in any case is simply a phe- 
nomcnological coupling constant in these calculations. In 
principle, the U that enters the coupling is different from 
the effective U in the denominator of x, Eq. (|29[) because 
of vertex corrections. Here, we will for simplicity ignore 
this distinction so we have just the basic RPA form of the 
spin-fluctuation interaction. A similar approach was used 
to discuss the microwave and thermal conductivity life- 
times |2lj . This approach will in principle include reso- 
nant spin excitations like those which have been proposed 
to explain the n resonance, i.e. the peak seen by neutron 
scattering at around 40 mcV in various cuprates, but we 
do not investigate these effects in detail here. Rather, we 
will examine what this approximation gives for the higher 
energy dependence of ri„ e i(u;, k| w=£k , T) numerically. 

In Fig. El we show the numerical evaluation of Eq. (|24p. 
as well as the RPA result with xo replaced by \ as 
in Eq. I|29[l . using the band parameters discussed above 
in Eqs. Q and (jSJ appropriate for an optimally doped 
cupratc. Here, the numerical integrations in Eq. (|24|l and 
(|2rJ)l are done using the technique described in Ref. [2(|. 
The value of U = 2. It chosen in the RPA case comes 
from early fits of spin-fluctuation theory to NMR data, 
and more recently to microwave and thermal conductiv- 
ity data on YBCO, but is expected to be reasonable for 
BSCCO as well because transport data suggest that the 
inelastic scattering rates at T c are very similar. The value 
of U = 6.7t in the case of the second order perturbation 
theory result was chosen so that the u> 3 behavior of both 
results matches for |u;|/Ao < 1 as shown in the inset. As 
noted, when \lu\ < A , the w-dependence of the inelastic 
scattering is determined by the phase space, and one gets 
an lo 3 dependence from a constant U interaction as well 
as from the spin fluctuation (RPA) interaction. How- 
ever, as shown in Fig. EI at higher energies the scattering 



reflects the frequency dependence of the effective inter- 
action. In the RPA case, the classical quasilincar energy 
dependence expected in spin fluctuation theory when the 
energy exceeds the spin-fluctuation energy is clearly vis- 
ible, and the crossover to the low-energ y a; 3 form takes 
place around 3Ao, as also found in Ref. [22 ]. 

In Figures El a) and b), we study the effect of temper- 
ature on the inelastic scattering. First, in Fig. El a) we 
consider the nodal scattering rate at w = 0, which col- 
lapses rapidly at T c due to the opening of the gap. In 
Fig. El b), the energy dependence of the scattering rate is 
plotted for various temperatures. It is interesting to note 
that for these band parameters, there is a considerable 
amount of upward curvature at low ui in r e i(o;, k| w=irk , T) 
even at T c , in contrast to the simpler nearly nested bands 
considered in polE^. 

IV. TOTAL SCATTERING RATE 

To include both types of scattering effects, we neglect 
interference processes between electron-electron collisions 
and impurity scattering entirely, and approximate the 
total scattering rate at the node by 

Ttot = r e l + Tinoi, (30) 

i.e. higher order processes like the influence of the inelas- 
tic scattering on the elastic scattering and vice versa are 
neglected. However, this is partially taken into account 
since we have chosen the parameters for the elastic and 
inelastic scattering on a phenomenological basis. 

In Fig. [7| we plot the on-shell total scattering rate r to t 
of a nodal quasiparticle as a function of both tempera- 
ture and energy. It is seen that the generic features, for 
reasonable assumptions about the magnitudes of the im- 
purity scattering rates and ranges, are as follows. As a 
function of temperature, one expects the nodal scattering 
rate to collapse at T c , and to obey a T 3 dependence at 
the lowest temperatures, whose coefficient is determined 
by both the elastic and inelastic processes. On the other 
hand, the low energy dependence of the scattering rate 
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FIG. 7: a) and b): Temperature dependence of the to- 
tal scattering rate r to t(0, k N , T)/A vs. T/T c for U = 2.2t 
and A /t = 0.2, and normal state elastic scattering rates 
of r o i(0, fcjv, T c )/Ao = 0.5 and 1, respectively. c) and 
d): Frequency dependence of the on-shell scattering rate 
r to t [lo, k\ UJ=£k )/Ao vs. uj/Aq for the same U and Ao and 
T = 0.1T C with r o i(0, fcjv,T c )/A = 0.5 and 1, respectively. 
Results for different values of K = 0.1 (solid), 1 (dashed), and 
K = oo (dotted) are shown. 



at low temperatures is dominated by the quasi-linear lo 
dependence of the elastic scattering. The generic result 
appears to be a quasilinear low-energy behavior for weak 
out of plane impurities, where the slope of this result is 
unrelated to the slope of the high-energy quasilinear be- 
havior determined by spin-fluctuation scattering. A very 
low-energy (lo ~ 1 — 2mcV) flattening or upturn due to 
unitary scatters G ~ lo) at low energies in Eq. ©, is 
expected to be present in principle. This contribution is 
likely to be negligible in BSCCO due to the small concen- 
tration of native planar defects in current samples, and 
the fact that the expected impurity bandwidth is com- 
parable to the energy resolution in ARPES experiments. 



V. CONCLUSIONS 

Understanding the lifetime of nodal quasiparticles is 
a fundamental problem of cupratc physics which must 
be solved if we are to agree that a d-wave BCS theory 
describes the optimally to over-doped superconducting 
state, the only part of the phase diagram apparently sus- 
ceptible to semiquantitative description at this time. In 
addition, the fact that nodal quasiparticles appear to be 
robust even for strongly underdoped samples may be an 
important clue to the physics of the pseudogap, and a 
theory of the nodal states within BCS which works at 
optimal doping may help us to understand this clue. Fi- 
nally, the lifetime of the nodal quasiparticles determines 



thermal and microwave conductivity as well as photoc- 
mission lineshapes and one would like to have a uni- 
fied description of these quantities within a single model. 
Here, we have calculated the on-shell lo — Sk nodal quasi- 
particle scattering rate which enters in determining the 
MDC lincwidth. Using a simple model which parame- 
terizes the forward elastic scattering in terms of a range 
k , a normal state scattering rate Tq, and an inelas- 
tic spin-fluctuation scattering parameterized by an RPA 
form with an effective Coulomb coupling U , we have stud- 
ied the temperature and w-dependence of the scattering 
rate. 

For a quasiparticle at the nodal momentum fcjv, as T 
decreases below T c , the elastic scattering rate decreases 
from its normal state value Tq. As the gap opens, it 
reaches a smaller value r e i(T) determined by the forward 
scattering parameter k and the renormalized gap velocity 
V2- This temperature dependence of the elastic scattering 
rate arises from the temperature dependence of the gap 
amplitude which controls the available phase space for 
scattering, suppressing it as the gap opens. Because the 
scattering rate depends exponentially on the gap, which 
opens rapidly in BSCCO, the clastic scattering rate has 
a cusp-like onset at T c , decreasing rapidly to its low tem- 
perature value. The inelastic contribution to scattering 
rate also decreases as the gap opens, with a cusp, and 
then decreases as (T/Aq) 3 at low temperatures. 

At low reduced temperatures, the w-dependence of the 
clastic scattering rate can exhibit a quasi-linear behav- 
ior, varying as ^ ("f\f ) if the scattering is forward or as 
To if the scattering is more isotropic. The inelastic 
scattering rate initially increases as w 3 /Aq so that there 
is an energy beyond which the inelastic scattering be- 
comes dominant. For energies greater than of order 3Ao, 
the inelastic scattering rate crosses over to a quasi-linear 
w-dependence with a slope of order one. 

Although the calculations we have presented here are 
straightforward, at the present time the experimental sit- 
uation regarding the direct measurement of the nodal 
scattering rate by ARPES is somewhat uncertain. There 
are some claims in the literature that the rate collapses 
in the SC state 0, as found theoretically here, and 
some that marginal Fermi liquid linear behavior consis- 
tent with a quantum critical point persists down to the 
lowest temperatures @ . As samples and resolution of the 
ARPES technique improve, we expect this discrepancy to 
be resolved and our prediction for the nodal quasiparticle 
MDC width to be testable. 



Acknowledgments 

The authors are grateful for discussions with 
D. Dessau, A. Fujimori, A. Ino, P. Johnson, and Z.- 
X. Shen. We would like to thank A. Chubukov for his 
comments regarding the second order on-shell inelastic 
scattering rate. Work was partially supported by ONR 



8 



grant N00014-04-0060 and NSF DMR02-11166. 



Y. Zhang, NR. Ong, P.W. Anderson, D.A. Bonn, 
R. Liang, and W.N. Hardy, Phys. Rev. Lett. 86, 890 
(2001). [14 
M. Chiao, R.W. Hill, C. Lupien, B. Popic, R. Gagnon, 
and L. Taillefer, Phys. Rev. Lett. 82, 2943 (1999). 
A. Hosseini, R. Harris, S. Kamal, P. Dosanjh, J. Preston, [15 
R. Liang, W.N. Hardy, and D.A. Bonn, Phys. Rev. B 60, 
1349 (1999). 

P.J. Hirschfeld, W.O. Putikka, and D.J. Scalapino, [16 
Phys. Rev. B 50, 10250 (1994). 

M.B. Walker and M.F. Smith, Phys. Rev. B 61, 11285 [17 
(2000). [18 
A.C. Durst and P.A. Lee, Phys. Rev. B 62, 1270 (2000). 

D. Dessau, J.D. Koralek, private communication; A. Fu- 
jimori, private communication; A. Ino, M. Taniguchi, pri- [19 
vate communication; P. Johnson, T. Valla, private com- 
munication. [20 
T. Valla, A.V. Fedorov, P.D. Johnson, Q. Li, G.D. Gu, 

and N. Koshizuka, Phys. Rev. Lett. 85, 828 (2000). [21 

A. Kaminski, J. Mesot, H. Fretwell, J.C. Campuzano, 

MR. Norman, M. Randeria, H. Ding, T. Sato, T. Taka- [22 

hashi, T. Mochiku, K. Kadowaki, and H. Hoechst 

Phys. Rev. Lett. 84, 1788 (2000). [23 

A. A. Kordyuk, S. V. Borisenko, A. Koitzsch, J. Fink, 

M. Knupfer, B. Biichner, H. Berger, G. Margaritondo, [24 

C. T. Lin, B. Keimer, S. Ono, and Yoichi Ando, 

Phys. Rev. Lett. 92, 257006 (2004). [25 

E. Abrahams and CM. Varma, Proc. Nat'l Acad. Sci. 
97, 5714 (2000). [26 
L. Zhu, P.J. Hirschfeld, and D.J. Scalapino, Phys. Rev. B 
70, 214503 (2004). 

H. Eisaki, N. Kaneko, D.L. Feng, A. Damascelli, 



P.K. Mang, K.M. Shen, Z.-X. Shen, and M. Greven, 
Phys. Rev. B 69, 064512 (2004). 

EW. Hudson, V. Madhavan, K. McElroy, J.E. Hoffman , 
K.M. Lang, H. Eisaki, S. Uchida, and J.C. Davis, Physica 
B 329, 1365 (2003). 

For a recent treatment of higher-order processes in the 
forward scattering li mit, see C.T . Rieck, K. Scharnberg, 
S. Schemer; cond-mat/0408320 

See, e.g. H. Won, S. Haas, D. Parker, S. Telang, A. Vany- 
olos, and K. Maki; cond-mat/0501463 
H.-Y. Kee, Phys. Rev. B 64, 012506 (2001). 
See also, G. Haran and A.D.S. Nagi, Phys. Rev. B 54, 
15463 (1996); ibid 58, 12441 (1998); M.L. Kulic and 
O.V. Dolgov, Phys. Rev. B 60, 
N. F. Berk and J. R. Schrieffer, 
(1966). 

D.J. Scalapino, 
1470 (1994). 

Hirschfeld, and D.J. 
224522 (2001). 
P.J. Hirschfeld, and D.J. 



S.M. Quinlan, 
Phys. Rev. B 49, 
D. Duffy, P.J. 
Phys. Rev. B 64, 
S.M. Quinlan, 



13062 (1999). 

Phys. Rev. Lett. 17, 433 



and N. Bulut, 



Scalapino, 



Scalapino, 



Phys. Rev. B 53, 8575 (1996). 

N. Bulut and D.J. Scalapino, Phys. Rev. B 45, 2371 
(1992). 

M.L. Titov, A.G. Yashenkin, and D.N. Aristov, 
Phys. Rev. B 52, 10626 (1995). 

J. Paaske and D.V. Khveshchenko, Physica C 341-348, 
265 (2000). 

T. Dahm and L. Tewordt, Phys. Rev. Lett. 74, 793 
(1995); Phys. Rev. B 52, 1297 (1995). 



